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Abstract 

For high energy processes (M 3> Aqcd) there are infrared safe hadronic shape vari- 
ables that have a calculable perturbative expansion in a s {M 2 ). However, nonperturbative 
power suppressed corrections to these variables are not well understood. We use the behav- 
ior of large orders of the perturbation expansion to gain insight into the nonperturbative 
corrections. Our results suggest that certain shape variables have nonperturbative correc- 
tions suppressed by fractional powers of Aq CD /M 2 . 
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Many large momentum transfer processes involving the strong interactions have a 
power series expansion in the strong coupling constant a s (M 2 ). In addition to this per- 
turbative expansion in a s , there are nonperturbative corrections suppressed by powers of 
Aq CD /M 2 , where M is the typical momentum transfer, M ^> Aqcd- I n some cases the 
nonperturbative corrections are well understood. For example, in R(M 2 ) = a(e + e~ — > 
hadrons) /a(e + e~ — > with M 2 = (p e - +p e +) 2 , the operator product expansion can 

be used to show that nonperturbative power corrections are proportional to the vacuum 
expectation values of gauge invariant local operators.0 The leading nonperturbative cor- 
rection is of order Aq CD /M 4 , and arises from the expectation value of the dimension four 
operator Tr G^C 1 " , where is the gluon field strength tensor. The operator product 
expansion shows that there are no corrections of order Aq CD /M 2 , since there is no gauge 
invariant dimension two operator that can be constructed out of quark and gluon fields. 
However, there are physical quantities that can be predicted using perturbative QCD for 
which the nonperturbative corrections are not well understood. For example, infrared safe 
event shape variables in e + e~ — ► hadrons have a perturbative expansion in a s (M 2 ). Since 
these variables weight the final hadronic states in a way that depends on their shape, they 
are not given by the imaginary part of the time ordered product of electromagnetic currents 
and an operator product expansion does not seem feasible. In this letter we examine some 
features of the nonperturbative corrections to hadronic shape variables in a toy model. 

The strong interaction coupling constant a s (M 2 ) falls logarithmically with the mo- 
mentum transfer M 2 . Nevertheless, perturbative QCD corrections can also produce power 
suppressed corrections to physical quantities, because the QCD perturbative expansion 
is an asymptotic series. Information on the power suppressed corrections to physical quan- 
tities can be obtained from an examination of the behavior of the perturbative expansion 
in a s at large orders. Consider a physical quantity P(a s ) that has a power series expansion 

oo 

P(a s ) = Y,Pna: +1 (1) 

n=0 

in the strong coupling a s (M 2 ). It's Borel transform is defined by 

B i p m = J2^r- ( 2 ) 

n=0 

1 There may also be non-perturbative instanton contributions to the coefficients of operators 

Q 
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The Borel transform B[P](t) may converge, even if the original series P(oe s ) is only an 
asymptotic expansion. This provides a definition of P(a s ) through the integral transform 

P oo 

P(a s ) = / dte-^ a ^B[P](t). (3) 



However, if the coefficients p n grow too fast with n, the Borel transform B[P] will have 
singularities that prevent using ([|) to define P(a s ). In an asymptotically free theory such 
as QCD, singularities in B[P](t) for positive t are associated with the infrared properties 
of Feynman diagrams, and are usually referred to as infrared renormalons. One can still 
obtain P(a s ) from B[P](t) using the inverse Borel transform eq. ([[]), provided one deforms 
the contour of integration in eq. (|3|) away from the real t axis, to avoid the renormalon 
singularity. The resultant expression for P(a s ) is no longer unique (or even real), and de- 
pends on whether the deformed path passes above or below the renormalon pole. Suppose, 
for example, that B[P]{t) has a simple pole at 

t = r , (4) 
where bo is the first term in the /3-function (negative in QCD), 

^=b al + O(al), (5) 

which governs the large M 2 behavior of the coupling constant in an asymptotically free 
theory, 

a s (M 2 ) -, (6) 

(-bo) In M 2 /A 2 



and u < is a constant (which is usually, but not always, an integer). The magnitude of 
the ambiguity in eq. (|3|) associated with this singularity can be estimated by the difference 
between the value of P obtained by deforming the t integration contour above and below 
the pole [Q. The difference is 

I API ~ e ~ u / b oa s 

< ^ (7) 
-(A^cd/M 2 ) ( - u) . 

Thus infrared renormalons produce power suppressed ambiguities in P. The infrared 
renormalon closest to the origin u = t = gives the dominant ambiguity. For P(M 2 ), 
the total cross-section for e + e~ — > hadrons defined previously, it is believed that the 
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perturbative expansion is not Borel summable, and has a renormalon at u = —2, which 
gives rise to an ambiguity in R that is of order Aq CD /M 4 m|. This ambiguity is of the 
same form as the leading non-perturbative correction in the operator product expansion 
for R(M 2 ), due to the vacuum expectation value of Tr G^ V G^ V . It is thought that there is 
an ambiguity in the definition of (0| Tr G^G^ |0), and that the renormalon ambiguity in 
the sum of the perturbation series cancels the ambiguity in the non-perturbative matrix 
element, to give a well-defined result for the total cross-section i?(M 2 )0. 

Renormalon ambiguities can only be absorbed into non-perturbative matrix elements 
if the contributions of both quantities to a physical observable such as R(M 2 ) have the 
same M dependence. Renormalon ambiguities can be absorbed in gauge invariant matrix 
elements provided the operator dimension is equal to — 2u. Explicit computations on 
renormalons are at a primitive level. Typically, one sums bubble graphs like those given 
in fig. |l| to determine the renormalon singularities @. Other graphs which have been 
neglected are presumably just as important, so the result of summing the bubble chain is 
not a proof of the existence of renormalon effects. Renormalon effects have been computed 
by Beneke || in a limiting case of QCD. One considers QCD with Nf flavors, and takes 
the limit Nf — > oo with a = Nfa s held fixed. Feynman diagrams are computed to leading 
order in a s , but to all orders in a. Terms in the bubble sum of fig. [1] with any number 
of bubbles are kept in this limit, since each additional fermion loop contributes a factor 
oc s Nf, which is not treated as small. The Nf — > oo limit has its limitations, since QCD 
is not an asymptotically free theory in this limit. The procedure used by Beneke is to 
write the Borel transform as a function of u = bot, (where bo is now positive), but still 
study renormalons for negative u. The singularities at negative u are then taken to be the 
infrared renormalons for asymptotically free QCD. This procedure was used by Beneke to 
study renormalons in R(M 2 ). It has been suggested by Brown and Yaffe || that there 
might be a renormalon at u = —1, which would lead to an order Aq CD /M 2 ambiguity 
in the perturbative expansion for R(M 2 ) that could not be compensated by the matrix 
element of a local operator. Beneke found by explicit computation that the renormalon 
singularity at u = — 1 vanished. 

In this letter we examine the behavior of shape variables at large orders in perturbation 
theory using a toy model. The model consists of a Z7i (1) x ^(l) gauge group (with gauge 
couplings e\ and e%)^ and a large number Nf of massless fermions / with charge (1,0). 
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We suppose there is a neutral scalar (j) of mass M that decays to the massless degrees of 
freedom through the dimension five interaction Lagrangian density. 



Ant = Ae ie2 F« (8) 



where A has mass dimension —1. We will compute the energy-energy correlation |T(| for 
<p decay to leading order in a\ and c*2, and to all orders in a = Nfa±. Physical situations 
analogous to this exist in QCD. For example a color singlet (Higgs) scalar (j) could decay 
to hadronic final states via the interaction £[ nt = Xg 2 (pTr G^G^. The non-abelian 
structure of QCD plays no role in the bubble chain sum of fig. [1], so we have simplified the 
computation by choosing an Abelian gauge theory, with two different £7(l)'s. 

The total (j) decay rate can be computed using the operator product expansion, and 
the computation is very similar to the computation of R{M 2 ). The leading operator is 
1, and the first correction comes from dimension four operators. One can also compute 
the Borel singularities using the methods of ref. H . It is straightforward to show that the 
total decay rate has no singularities in the Borel plane. 

A more interesting computation is that of event shapes in <p decay. There is no 
operator product expansion for the event shapes, so it is interesting to see whether there 
exist renormalon singularities, and whether they are at integer values of u. A convenient 



characterization of event shapes in <fi decay is the energy-energy correlation |T(J 
V(cosx) = 

El -n=, — 7^=, — ; 7T~ ( tC) ( \ 5 (cos 9a — cosy) dEidE^d cos 6>i- 
J dE t dEj dcosOij \M J \M J K 3 

+ , (1 _ rosx)E |j_(!) 2 d£( . 



(9) 



In eq. (^), the double sum is over all pairs ij of particles in the final state with i ^ j, 
and ij and ji are both included in the sum. The angle between the three-momentum 
vectors of particles i and j in the <fi rest frame is denoted by 0^, and Ei represents the 
energy of particle i. The Fox- Wolfram moments [JTTJ Hl are obtained as integrals of the 
energy-energy correlations with Legendre polynomials, 

l 

dcos x V (cos x) Pl (cosx) = Y Hl, (10) 

l 



where V is the total <fi decay rate. Integrals of P(cosx) against any "smooth" weighting 
function w(cosx) are infrared safe event shape variables. 

The Feynman diagrams of fig. contribute to the (p decay distribution to lowest order 
in A, <X\ and az, and to all orders in a = NfOt,\. We will restrict our analysis to shape 
variables that can be obtained by integrating V against a weighting function that vanishes 
at cosx = ±1. This simplifies the calculation because the Feynman diagram in fig. (2a) 
and the single sum in eq. (^) don't contribute. The remaining Feynman graphs are easy 
to compute, and one obtains 



P(cos X )=/ dq 2 ^^, (11) 



2 dV(cosx) 

r icos x) = / «C 
Jo 

where 

= 2mx2 ^ i F { f COS v) (12) 

dq 2 N f |l + n(^)|^ [q ' SXj ' { } 

Here q 2 is the invariant mass of the fermion pair, 

U(q 2 ) = -a b In (-q 2 /A 2 ) , b = 1/3tt, (13) 

is the fermion bubble contribution to the gauge propagator, and 

q 2 (l — g 2 ) 4 (1 + cos x) 2 + <? 4 (1 — cos xY 



F (q , cos x) 



+ 



[(1 + cosx) +q 2 (1-cosx)] 
9 ((1 + q 2 ) 2 (1-cosx) -Sq 2 ^jq 2 [(1-cosx) (l + <f) - 4g 2 ] 

7 /2 2 "I ^ 

(1- cosx) (l + <7 2 ) (l-cosx)-8g 



(14) 



with g 2 = q 2 /M 2 . In eq. ([14]) terms proportional to 5(1 + cosx) an d 5(1 — cosx) are 
neglected. The quantity A in eq. (|T3| ) is proportional to the subtraction point \i used to 
define the finite part of the vacuum polarization, IT, and is in general subtraction scheme 
dependent. 

We are interested in the Borel transform of V (cos x) with respect to the variable 
a = NfOt\. The dependence of V (cosx) on a is i n the function 



9(a) 



a 2 a 2 



,1 + U(q 2 )\ 2 (1 - ab In (q 2 /A 2 )) 2 + b 2 a 2 n 2 ' 

StStW E + 1) 2m + 1)6^ 2 (In (, 2 /A 2 ))— . 

( ([15]) (|T3|) equation 1 hequation 1515) 



(2m + 1 

m=0 n=2m 



The Borel transform of g (a) is 



B fa] w = E t^ttW E -k +2ra <" +2m+1 (m ( 9 2 /a 2 ))" , 

m=0 ^ n=0 

= ( ? 2 /A 2 ) M sm(vr6 t), 
= (q 2 /A 2 ) u sin (ttm) , 

((|15l)(ll^egtta£zorel6egm£zcml616) 

using u = bot. 

Consider an infrared safe shape variable S v that corresponds to a weighting function 
w(cosx) that behaves like w(cosx) ~ (1 + cosx) v near cosx = —1. The singularities 
in the Borel transform of this shape variable that arise from the integration region near 
cosx = — 1 are determined by the value of F(q 2 , cosx) for small q 2 and cosx near —1. In 
this region 

q 2 (1 + cosx) 2 + 4g 4 
F(g 2 ,cosx) ~ — ; — 5^-- ( (p!5[) (|T5|) equationl 7 equation 1717) 



[(1 + cosx) +2<? 2 ] 

Using eqs. (]TT|) — we find that the Borel transform of has infrared renormalons. For 
rj < 1 the singularity closest to the origin is a simple pole at 



u 



-n =^> t = —r]/bo. (tfnfy($nfyequationl8equationl818) 



For r\ = 1 there is potentially a simple pole at u = — 1. However the sin(7i"u) factor cancels 
the pole at u = — 1, and there is no singularity in -BfS'i] at u = — 1. The singularity closest 
to the origin is a simple pole at u = —2. 
An integral of the form 

dq 2 h (q 2 ) B [g] = sin (iru) / dq 2 h (q 2 ) (q 2 /A 2 ) " 

( (HH) (HI) e 1 ua ti on 1 9eg uat i on 1 9 1 9 ) 
determines the singularities in B[S V ]. If h(q 2 ) ~ (q 2 ) k for small q 2 , there is no singularity 
provided k is an integer, and one obtains a simple pole singularity at u = — k— 1 if k is not an 
integer. There is a simple pole singularity at u = —k — 1 for integer k if h(q 2 ) ~ (Q 2 ) k m (q 2 ) 
for small q. The weighting function (1 + cosx) v produces an effective (q 2 ) V 1 behavior! 
for h(q 2 ) at small q 2 , and hence has a simple pole singularity at u = —rj. For n = 1 (or 



q uation20footnote22 Although h ( q 2) is s i ngu l a r as q 2 -» 0, it is integrable so that 5, is infrared safe. 
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more precisely, for w (cosx) = cos 2 x — 1 which vanishes at cosx = ±1), h(q 2 ) in eq. ( fl^ ) 
has the form cq + c\ g 2 lng 2 for small q 2 , where cq and c\ are constants. The absence of 
a lng 2 term (which is infrared safe) implies that there is no renormalon at u = —1. The 
q 2 \rtq 2 term produces the pole at u = —2. In more realistic examples, one expects that 
there will be a lng 2 term, and hence a renormalon singularity at u = — 1 as well. 

The behavior of QCD hadronic event shape variables can be inferred from the toy 
example using the method of Beneke ||. Our results imply that Fox- Wolfram moments 
will have renormalons at integer values of u. Event shape variables defined by weighting 
functions that are polynomials in cos x (e.g. Fox- Wolfram moments) have Borel transforms 
that are more singular than the Borel transform of the total rate. More interesting are the 
shape variables S v (0 < r\ < 1) which have a renormalon at u = —rj. These variables have 
ambiguities of order (Ag^/M 2 )^ from the high-order sum of the perturbation expansion. 
Non-perturbative power suppressed corrections are expected to be of the same order and 
cancel this ambiguity, leaving a residual power suppressed correction. The unusual power 
dependence of the shape variables S v arises from weighting functions w(cosx) that are 
continuous but not differentiable at cosx = — 1- O ur results indicate that hadronic shape 
variables which are infrared safe but do not treat the region where infrared divergences 
cancel between Feynman diagrams in a very smooth fashion can have large power sup- 
pressed corrections. We hope to present results for more realistic situations in a future 
publication. 
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Figure Captions 

>Dequation20figurell. Graphs that give an infrared renormalon in QCD. 

)|)equation20figure22. Feynman diagrams contributing to the event shape distribution to lowest order 

in A, «i and «2, and to all orders in Nfa±. 
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Figure 1 



Figure 2 



